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Abstract 

We describe the universe as a local, inhomogeneous spherical bubble embedded in a flat 
matter dominated FLRW universe. Generalized exact Friedmann equations describe the ex- 
pansion of the universe and an early universe inflationary de Sitter solution is obtained. A 
non-perturbative expression for the deceleration parameter q is derived that can possibly de- 
scribe the acceleration of the universe without dark energy, due to the effects associated with 
very long wave length super-horizon inflationary perturbations. The suggestion by Kolbe et 
al. 9 that long wave length super-horizon inflationary modes can affect a local observable 
through inhomogeneities is considered in the light of our exact inhomogeneous model. 

e-mail: j moffat @perimeterinstitute . ca 

1 Introduction 

In a recent article, we investigated a cosmology in which a spherically symmetric perturbation 
enhancement is embedded in an asymptotic FLRW universe pQ. The perturbation enhancement 
is described by an exact inhomogeneous solution of Einstein's field equations. We found that 
the large-scale inhomogeneities can lead to a reinterpretation of the luminosity distance of a 
cosmological source in terms of its red shift z, owing to the observer dependence of these quantities. 
The time evolution and the expansion rate of the inhomogeneous universe can lead to intrinsic 
effects such as cosmic variance at large angles and long-wavelength perturbations not described 
by a FLRW homogeneous and isotropic universe. Therefore, the interpretation of the data using 
a FLRW model that the accelerating expansion of the universe is caused by dark energy may be 
misleading. This is important, for it is difficult to explain theoretically the postulated dark energy 
that causes the acceleration of the universe. The model also leads to an axis pointing towards the 
center of the spherically symmetric large scale perturbation enhancement with dipole, quadrupole 
and octopole moments aligned with the axis. It was shown that the luminosity distances and red 
shifts observed by different observers located at spatially different points of causally disconnected 
parts of the universe can have varying values. A spatial average of all these observations leads 
to an intrinsic cosmic variance in e.g. the deceleration parameter q. The distribution of CMB 
temperature fluctuations can be unevenly distributed in the northern and southern hemispheres. 

The popular explanation for the observed large-scale homogeneity of the universe is that the 
universe underwent an initial inflationary period with more than 60 e- folds |2J. The inflationary 
cosmic expansion can stretch an initially small, smooth spatial region to a size larger than the 
horizon size today and explain the present day large-scale homogeneity. The question arises as 
to whether the initial local patch can be sufficiently homogenized to allow inflation to begin 
In the following, we shall consider the universe as an expanding bubble with an inhomogeneous 
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metric and generalized Friedmann equations, including inhomogeneous density and pressure and 
a cosmological constant. Our main assumptions arc spherical symmetry and an inhomogeneous 
barytropic fluid that satisfies an equation of state. For the case of a spatially fiat inhomogeneous 
early universe, we obtain a de Sitter inflationary solution. 

The acceleration of the expansion of the universe deduced from Type la supernovae observations 
and the CMB WMAP data |SJ EI has been interpreted as due to the cosmological constant 
(vacuum energy), modifications of Einstein's gravitational field equations at large distances |7], 
and quintessence fields 8 . The quintessence explanations postulate a new form of matter with 
negative pressure called dark energy. Recently, it has been suggested that the acceleration is caused 
by very long wavelength, super-horizon perturbations generated by a period of inflation in the early 
universe |5J EH . The backreaction of perturbations on an FLRW background universe has been 
the subject of investigation by several authors The predictions based on perturbation theory 
are limited by the condition <E> 1, where $ is the gravitational potential. 

The inflationary perturbation modes whose wavelengths presently are smaller than A < 10 
Mpc have entered the non-linear regime and have generated galaxies and clusters of galaxies, while 
longer wavelength modes at super-horizon scales > c/ H are entering the linear regime today. 
The effects of the sub-horizon modes are small due to fact that Sp/p ~ 10~ 5 at the surface 
of last scattering. Therefore, these sub-horizon modes produce negligible corrections at second 
order ~ 10 -8 . However, the super-horizon modes could potentially create a correction to the 
deceleration parameter q, large enough to remove the need for dark energy. It has been argued 
recently |12M13llT%lll5| that second order perturbation effects of the form <(>V 2 $ are described by 
a renormalization of the local spatial curvature and cannot (for a positive energy density) produce 
a negative deceleration parameter. 

One problem with the perturbation calculations is that they ignore all higher gradient terms 
V"<I>, and any non-perturbative effects that can have a significant influence on the inhomogeneity 
contributions due to very long wave length modes at super-horizon. These effects will occur for 
inflationary models in which the number of e-folds of inflation is much larger than the 60 e-folds 
required to create physically satisfactory fluctuations. In the following, we shall use the exact 
inhomogeneous model of ref. PP to derive a formula for the deceleration parameter q that is non- 
perturbative and whose variance can lead to a negative value for q without dark energy and a 
cosmological constant A. 



2 Inhomogeneous Friedmann Equations 

Our action takes the form 

S = S G + S M , (1) 

where 

^ = I^/<*W^(*-2A). (2) 



The matter action is given by 



s ,,= I d^V^gi^d^d^ - v (</>)], (3) 



where <f> is a scalar matter field and V(<j>) is a potential. 

For the sake of notational clarity, we write the FLRW line element 

dr 2 



ds 2 = dt 2 - a 2 (t) I - + rWj , (4) 

where d£l 2 — d9 2 + smd 2 d(j> 2 . The general, spherically symmetric inhomogeneous line element is 
given by [JJ1 UH DB [U EDI EJ HI HJ : 

ds 2 = dt 2 - X 2 (r, t)dr 2 - R 2 (r, t)dn 2 . (5) 
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The energy-momentum tensor T M „ takes the barytropic form 

n = (p+pK^-p^, (6) 

where u M = dx^/ds and, in general, the density p = p(r, t) and the pressure p = p(r, t) depend on 
both r and t. We have for comoving coordinates u° — 1, u % = 0, (i = 1, 2, 3) and g liv u il u v = 1. 
The Einstein gravitational equations are 

Gfj, u = Rp, u - -g^vll + Ag^ = -SnGT^, (7) 

where 1Z = g^R^ and A is the cosmological constant. Solving the Goi = equation for the 
metric ©, we find that 

*(-«>- w 

where i?' = dR/dr and /(r) is an arbitrary function of r. 
We obtain the two generalized Friedmann equations pQ: 

i? 1 R 2 11, , 2 , lR'R lR"ff 4irG, „ . 1. 
i? + 3i^ + 3i^ (1 - ^ - 3 + 3 W = " — ( " + 3P) + 3 A ' (10) 
where i? = dR/dt. 

3 De Sitter Inflationary Solution 

Let us now consider the very early universe and retain the pressure p and the cosmological constant 
A. We shall picture the early universe as an expanding spherically symmetric bubble with its origin 
at the big bang. We choose f(r) — 1 for all values of r, and obtain the generalized Friedmann 
equations 

| +2 || =8lGp + A, (11) 
R 1 R 2 1R' R 4ttG . „ , 1 . 

R + 3^-3R^R = - — (P + *P)+3 A - ( 12 ) 
By using the notation H± = R/R and H r = R! /R', we can write {TJ and 112|) as 

H\ + 2H x H r = 8nGp + A, (13) 

| + \Hl - l -H r H^ = -*£(p + 3p) + 1a. (14) 

For H±(r,t) = H r (r,t) = H(t) — a(t)/a(t) and R(r,t) = a(t), we obtain the Friedmann equations 
of FLRW for a spatially flat universe: 

H 2 = ^ + \k. (15) 

- = ~{p + 3p) + U (16) 
a 3 3 

where p = p(t) and p = 
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From the Bianchi identities V ' V G^ V — 0, we obtain for the inhomogeneous model the conserva- 
tion law 

V„T^ = 0. 



This becomes 



1 



-9 



-.d v [V ZI 9{p+p)u ll u u ]+Yi x v xu v u x = 0. 



For T^oo = and y/^g = X(r, t)R 2 (r, t) sin#, we obtain 

dp Id 



dt XR 2 dt 



[XR 2 (p + p)}. 



From JBJ) we get 



P* = \P-\$<I>) 2 -V{<I>). 
Here, <fi = (j>(r, t) depends on both r and t. From l|18|) we obtain 

1 d 



™ ' XR 2 dt 
Differentiating l|2U|) with respect to t we have 

•■ 1 d 



{XR 2 ){p 4> +p dp ) = Q. 



-<W + 4iV(*y d/dt{XR2) 



Let us assume an equation of state 



XR 2 



p(r,t) = wp(r,t), 



r (V0) 2 /</>] = o. 



(17) 
(18) 

(19) 

(20) 
(21) 

(22) 

(23) 
(24) 



where w is a constant. Then, a vacuum equation of state has w = — 1, while the matter dominated 
universe has w = 0. We choose for simplicity /(r) = 1 for all values of r and we assume that the uni- 
verse is dominated by the cosmological constant with p(r, t) = —p(r, t) and \/A/3 = ^/87rGp vac /3. 
For the special case of a spatially flat universe, a solution to Eas. l|ll|l and (|12|) is given by 



R(r, t) = R exp 



By a change of radial coordinate 



Rq exp 



(v^V3)(r + t) 



(0V/3)r 



we obtain the isotropic and homogeneous de Sitter metric: 



ds = dt — exp 



2(VA/3)* 



=2 



exp 



2(^/173) 



(25) 



(26) 



(27) 



Thus, the inhomogeneous expanding bubble can inflate when the cosmological constant (vacuum 
energy) dominates the expanding bubble. The solution 125JI is a special solution and a more 
general inhomogeneous inflating solution should exist for f(r) ^ 0, but it does illustrate that an 
inflationary solution of Einstein's field equations exists for our more general inhomogeneous metric. 

Assuming that the matter scalar field <f> dominates in the early universe, and substituting for 
p,p and Pcf, in (f 1 If) and l|12|l . we obtain for A = 0: 



R2 n RR ' a „ 

R 2+2 RR'=^ G 



^2 + l ( ^ + F(0) 



(28) 
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R 1 R 2 1 R' Ft 8itG . ; 5 „, ... /„„> 
5 + 3^-3^5 = — (29) 
We can obtain an inflationary solution of the form H25|) by assuming that V(<f>) ^> 2 , (V) 2 ~ 
and V((f>) ~ constant. However, we are required to make assumptions about the scalar field (f>, 
namely, that in a local patch of inflation it only depends on time t, so that the wave equation (|23[l 
can be solved to give the sub-horizon quantum fluctuations that seed the growth of structure. This 
appears to be an apparently unavoidable fine-tuning of the primordial scalar matter fields that 
occurs in generic inflationary models. 



4 Late-Time Matter Dominated Universe 

The late-time matter dominated universe will be pictured as a large scale perturbation enhance- 
ment that is described by an exact inhomogeneous spherically symmetric solution of Einstein's 
field equations. The perturbation enhancement is embedded in a matter dominated universe that 
approaches asymptotically an Einstein-de Sitter universe as t — > oo. An observer will be off-center 
from the origin of coordinates of the spherical perturbation enhancement. 

For the matter dominated Lemaitre-Tolman-Bondi (LTB) [1611171 ITS) model with zero pressure 
p = and zero cosmological constant A — 0, the Einstein field equations demand that R{r,t) 
satisfies 

2RR 2 + 2R{1 - f 2 ) = F(r), (30) 

with F being an arbitrary function of class C 2 . There exist three possible solutions depending on 
whether / 2 < 1,= 1, > 1 and they correspond to elliptic (closed), parabolic (flat), and hyperbolic 
(open) cases, respectively. 

The proper density of matter can be expressed as 

F' 

9 = WttGR'R 2 ■ (31) 

We can solve l|3*T)) to obtain 

p c FF 2 B \ R 2 RR' ) y ' 

where 

H 2 S = H\ + 2H ± H r , (33) 
is an effective Hubble parameter and we have 

^ Pc =- 2+ 2--. (34) 

We have the three possibilities for the curvature of spacetime: 1) f 2 > 1 open (SI — 1 < 0), 2) 
f 2 = 1 flat (n - 1 = 0), f 2 < 1 closed (SI — 1 > 0). 

Since the WMAP data 6 shows that the universe is spatially flat to within a few percent, we 
shall consider the globally flat case f 2 (r) = 1. The metric reduces to 

ds 2 = dt 2 - R 2 (r, t)dr 2 ~ R 2 (r, t)dQ 2 . (35) 
A solution for the matter dominated (w=0) universe is 

R(r,t) = r{t + P(r)} 2 /\ (36) 
and f3(r) is an arbitrary function of r of class C 2 J5j . The metric l)35|l becomes 

ds 2 = dt 2 - (t + l3) 4 / 3 (Y 2 dr 2 + rW), (37) 
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where 



and 



Y = 1 



2r/3' 
3(t + (3) ' 



1 



(38) 



(39) 



P 6irG{t + fj) 2 Y' 

The arbitrary function can be specified in terms of a density on some spacelike hypersurface 
t = U). The metric and density are singular on the two hypersurfaces t + /3 = and Y = 0, 
namely, t\ = —ft and £2 = —0 — 2r/3'/3, respectively. The model is only valid for t > S(r) = 
Max[£i(7"),t2(r)], and the hypersurface t(r) = E(r) defines the big-bang. However, our pressureless 
model requires that the surface t(r) = S(r) describes the surface on which the universe becomes 
matter dominated (in the LFRW model this occurs at z ~ 10 4 ). We observe that even in the 
spatially flat LTB model, different parts of the universe can enter the matter dominated era at 
different times. Our expanding inhomogeneous, spherically symmetric bubble is embedded in the 
matter dominated metric (|37|l . For f3 = and in the limit t — > 00 we obtain the Einstein-de Sitter 
universe 



ds 2 = dt 1 



2 (t)(dr 



l dn 2 ), 



(40) 



where a(t) — t 2 / 3 . Thus, for [3 — we obtain the FLRW model. Moreover, the expanding flat 
LTB model necessarily evolves to the homogeneous and isotropic FLRW model for a non- vanishing 
density, whatever the initial conditions. 



5 The Inhomogeneous Cosmology Deceleration Parameter 

Let us expand R(r, t) in a Taylor series 



R{r, t) = R[r, t - (t - *)] = R(r, to) 



R(r, t ) 
R{r,t ) 



1 - (to - t)H ± -^0- t) 2 q(r,t Q )H 2 ± 



= R(r, to) 

where to denotes the present epoch and H ± = R(r, to)/R(r, to). Moreover, we have 

R(r,t )R(r,t ) 



q{r,t ) = - '- 



By substituting for R from Eg. 1121) . we obtain 



q(r, to) 



47rp (r, t ) 



R 2 (r,t ) 



A 



1 H 0r (r,t ) 



3 3H 2 ± (r,to) 3H 2 ± (r 7 t ) 3 H 0± (r,t )' 



(41) 



(42) 



(43) 



where H Qr {r,t Q ) = R'(r,t )/R'(r,t ). 

If we set Ho±(r,to) — Ho r (r,to) = H(to) where H(to) = a(to)/a(to), then we obtain the 
spatially flat FLRW expression for the deceleration parameter: 



3H 2 



n — ^ A ■ 

3H 2 2 A 



(44) 



We see from (|43H that different observers located in different causally disconnected parts of the 
sky will observe different values for the deceleration parameter q, depending upon their location 
and distance from the center of the spherically symmetric perturbation enhancement. This can 
lead to one form of cosmic variance, because the spatial average of all the observed values of local 
physical quantities, including the deceleration parameter q, will have an intrinsic uncertainty. 
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Let us rewrite (14311 in the form 



q = g +fiom± - ^OA-L - &0H.L, (45) 

where 

Anpo(r,t ) A 1 H 0r {r,t ) 

3 ^o±( r ^o) 3JJ^ x (r,t ) 3 J ff 0X (r,i ) 

The variance of q is given by the exact non-perturbative expression: 

v^( q )^(q 2 -(q)Y^ = q, (47) 

where 



denotes the ensemble average. 

We now set A = in (|45|) and obtain 



q = g +«0mi - Ooh±- (49) 

If we have 

^Offl > g + ^Oml, (50) 

then the cosmic variance var(<7) can be negative and cause the universe to accelerate without 
a cosmological constant or dark energy. If the super-horizon long wave inflationary modes can 
generate non-perturbatively a sufficiently large effect of O(l), then this can drive var(g) to negative 
values. 



6 Conditions for Deceleration Parameter 

Hirata and Seljak ^3] have recently raised some critical questions regarding the Kolbe et al. [HI 
suggestion that long wave length super-horizon perturbations can cause q to be negative and not 
violate the strong energy condition p + 3p > 0. They used non-FLRW methods to discuss the 
deceleration parameter q. They begin by defining a local Hubble value 



Hi = -V^tt", (51) 



and a deceleration parameter 



9i = -1 - -^ V m#i- (52) 

Here, u M V M is the Lagrangian proper time derivative d/t plopeT associated with a matter particle. 
The expansion tensor M „ is related to H 1 by 9 = 9^^ = 3i?x- The matter particles follow geodesies 
with O^u" = 9^ v u^ = 0. From the Raychaudhuri equation 

d9 9 2 

-j t = - j - + - R^vPu", (53) 

where R^ v is the Ricci tensor and a and to denote the shear and vorticity, respectively, it follows 
that 

Hhi = ~ vWn) - \R^u v . (54) 
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From Einstein's field equation J7J for A = we obtain 

Hfqi = \{o» v a% - ufvufv) + ^-(p + 3p). (55) 

We have ' v cr v ' M > 0, so that if the strong energy condition is satisfied, p + 3p > 0, and the 
vorticity lo^u = 0, it follows that q\ > 0. Hirata and Seljak now state that in the synchronous 
comoving gauge used by Kolbe et al. to calculate the perturbations up to second order, u) = 0, 
and the deceleration parameter must satisfy q\ > when the strong energy condition is satisfied. 

In our model the vorticity ui ^ for our exact inhomogcncous solution that describes the large 
scale perturbation enhancement, so it does not necessarily follow that qi > 0. We do not have 
to use the synchronous gauge to do calculations of the non-pertubativc inhomogcncity and so, in 
general, u> ^ 0. 

Harati and Seljak show that for a single inflaton inflationary model, the long wave length 
fluctuations at super-horizon cannot give rise to an effect O(l) for var(g). However, this result is 
model dependent and a more general hybrid inflation model could lead to infrared modes that can 
lead to a shifting of var(g) to negative values, and explain the acceleration of the universe without 
negative pressure dark energy or a cosmological constant. 

If the vorticity uj ^ and the shear a ~ 0, then we can consider the possibility that the long wave 
length super-horizon inflationary modes can lead to a negative q. However, it can be argued on the 
basis of standard inflationary models that due to the large increase of entropy during re-heating, 
there will be a large increase of p + 3p in (|55|l . while the vorticity is not increased accordingly 
and, therefore, one may not expect that q can be driven to negative values by a sufficiently large 
vorticity value. But this depends on the re-heating model and may not be a generic feature of 
inflation. 

Alternatively, one could investigate these issues in a non-inflationary model, such as the variable 
speed of light bimetric gravity theory [231 , which does predict a scale invariant spectrum with the 
spectral index n s ~ 0.97 in agreement with the WMAP result. The superluminal mechanism 
that solves the horizon and flatness problems, namely, that the ratio 7 = c 7 /c g , where c 7 and c g 
denote the speed of light and the speed of gravitational waves, respectively, becomes unity soon 
after it reaches a high value in the very early universe, as the effect of a phase transition with 
c g = c = constant (c the currently measured speed of light). The fluctuations that form the seeds 
of large scale structure are born super-horizon but they do not lead to the same consequences as 
inflation with a large number of e- folds. 

7 Conclusions 

We have shown that a special de Sitter inflationary solution exists for a spherically symmetric 
and inhomogeneous expanding bubble. Further work must be carried out to find more general 
inhomogeneous de Sitter- like solutions from our generalized Friedmann equations. 

We have argued that an exact non-perturbative treatment of inhomogeneities in a cosmology 
describing a spherically symmetric expanding universe, embedded in an asymptotic Einstein-de 
Sitter universe, can significantly modify local observables such as the red shift, expansion rate and 
the luminosity distance when super-horizon, long wave length modes are taken into account. We 
expect that this will hold true even in the case in which only adiabatic modes are present and 
spatial gradients and non-perturbative effects are not ignored. 

Care must be taken when arguments are made about the sizes of effects associated with the 
super- horizon inflationary modes in cosmology. According to our results, the various features of the 
power spectrum derived from the WMAP data could also be subject to significant modifications at 
the non-perturbative level depending on the model of inflation adopted. The effects of an expansion 
parameter 6 and shear and vorticity parameters should be reinterpreted in our inhomogeneous 
model together with the effects of very long wave length super-horizon inflationary modes, before 
conclusions are drawn about their relative importance when compared to an FLRW universe. 
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The fact that we should consider the effects of super-inflationary modes in cosmology and that 
they lead to a cosmic variance for various local physical quantities, even though these modes cannot 
presently be observed, may be inescapable and lead to an intrinsic uncertainty in our interpretation 
of the universe. 
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